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We analyze the prospects for stabilizing Majorana zero modes in semiconductor nanowires that are proximity- 
coupled to higher- temperature superconductors. We begin with the case of iron pnictides which, though they 
are s-wave superconductors, are believed to have superconducting gaps that change sign. We then consider the 
case of cuprate superconductors. We show that a nanowire on a step-like surface, especially in an orthorhombic 
material such as YBCO, can support Majorana zero modes at an elevated temperature. 
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I. INTRODUCTION 

Majorana zero modes have been predicted to occur in cer- 
tain fractional quantum Hall states 1 5 and in superconduc- 
tors in which time-reversal and spin-rotational symmetries are 
broken 6 19 . Very recently, exciting experimental progress^^l 
has been made with semiconductor nanowire d 16 ! 17 * in which 
superconductivity has been induced by proximity to an s-wave 
superconductor (see, in addition, interesting results on a topo- 
logical insulator-superconductor interface in Ref. |24| |. The 
single-particle gap in the nanowire, which is the energy scale 
that controls the stability of the zero mode, is determined by 
the spin-orbit coupling in the wire, the applied Zeeman field in 
the direction of the wire, and the superconducting gap. While 
the first of these three quantities is determined primarily by 
the nanowire and the interface, second and third are limited 
primarily by the superconductor. Therefore, it is natural to 
consider superconductors with higher- T c s. In this paper, we 
focus on the pnictide and cuprate superconductors. 

The principal complication with unconventional supercon- 
ductors such as the pnictides and cuprates, vis-a-vis the prox- 
imity effect, is that the superconducting gap changes sign 
from one part of the Brillouin zone to another. Therefore, 
it is not even clear that such a superconductor can induce 
a non-zero superconducting gap in a nanowire. In addition, 
the hole-doped cuprate superconductors, at least, have gapless 
bulk excitations, which would hybridize with a Majorana zero 
mode, thereby giving it a finite lifetime. However, if these 
issues could be circumvented, then a higher- T c superconduc- 
tor could lead to a Majorana zero mode that is protected by a 
larger gap, A. This, in turn, can lead to a substantially smaller 
splitting e ~ LA / Vp or decay rate e~ A / T for zero modes. 

In this paper, we show that it is possible for a bulk s± su- 
perconductor (such as, according to many theories, a pnictide) 
to induce a large superconducting gap in a nanowire. As a re- 
sult of this gap, a Majorana zero mode is stabilized at each 
end of the wire. We further show that a d-wave supercon- 
ductor can induce a large superconducting gap in a nanowire 
that is coupled to it on a step-like terraced surface. If the 
superconductor has d + id pairing symmetry, then the wire 
would support a Majorana fermion zero mode at each end. 
(There is no known superconductor with such a gap symme- 
try, although one might expect, on symmetry grounds, that an 



id xy component would be induced in a d x 2_ y 2 superconduc- 
tor by the application of a magnetic field.) If the superconduc- 
tor has d x 2_ y 2 pairing symmetry, then the Majorana fermion 
zero mode may decay into the bulk. We compute this decay 
rate and analyze the conditions under which the zero mode is 
stabilized. 

For a recent discussion from a different perspective of 
topological superconductivity in semiconductor-cuprate struc- 
tures, see Ref. |25j For interesting recent experimental results 
finding a large superconducting gap induced in a topological 
insulator by proximity to a cuprate superconductor (as pro- 
posed in Ref. |26l . see Ref. |27] 



II. BASIC SETUP 
A. Action 



We assume that our system consists of a quasi- ID semicon- 
ductor nanowire of length L x aligned along the x-direction, 
with width w in the y-direction and width w z in the z- 
direction. We will assume that w z <C w <C L x and ne- 
glect motion in the z-direction but allow for the possibil- 
ity of multiple sub-bands in the y-direction. We will con- 
sider two nanowire widths: w — 50 nm, which is techno- 
logically possible (e.g. by e-beam lithography), as well as 
w = 100 nm, which is the width in recent experiments^^ 
on nanowires coupled to conventional s-wave superconduc- 
tors. The nanowire is described by the action Smw while 
the superconductor is described by the action Ssc\ these are 
shown below. The nanowire is lying on a (001) surface of 
the superconductor, and interactions between them are de- 
scribed by St- We will assume that the nanowire is pointing 
in the (100) direction. In the case of the cuprates, this will 
be the optimal arrangement (as will a nanowire pointing in 
the (010) direction, which is equivalent). If the wire points in 
the (1 10) direction, then it will be aligned with the nodes and 
the induced gap will generically vanish. The effective action 
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S = Snw + Ssc + St can be written in the form: 
2vr 
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Here, the index m is a Matsubara index, and i is a sub- 
band index in the nanowire. In Eq. ([TJ, we have suppressed 
the spin and particle-hole indices on W and c. We define 

^m,i,k x = (V'm,i,fc x 'f , "0m,i,fc m ^i V 1 — m,i. — fc x ^> V 1 — m,i, — fe^f) ! 

where ipm^.k^a is the annihilation operator for an electron of 
Matsubara frequency cj m , momentum /c x , and spin a in the 
i th sub-band of the nanowire. c m q is defined analogously as 
a spinor in spin and particle-hole space. 

The functions ti(k x ,q) are the tunneling amplitudes be- 
tween the superconductor and the i th sub-band of the semi- 
conductor nanowire. We assume for simplicity that the elec- 
tron wavefunction in the i th sub-band of the nanowire takes 
the form ip(x,y) = y^sin(^) tp(x). (The normaliza- 
tion is for later convenience.) The kinetic energy of an iso- 

h 2 k 2 

lated nanowire is e^ T = 2m^' anc ^ m * 7 ^ x ~ 9^bB x /2, 
(i, and a are, respectively, the effective mass, Zeeman split- 
ting, chemical potential, and strength of spin-orbit coupling 
in the nanowire. For an InSb nanowire, typical values are 
to* = 0.015 TO e , g « 50, and a = 200 meV A. We will 
assume that the chemical potential [i = several meV, and is 
tunable with gate voltage. In the absence of superconductiv- 
ity, the chemical potential is related to the Fermi momentum 
k,F,i in the i th sub-band according to 
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In this paper, we consider superconductors with higher T c : 
cuprate and pnictide superconductors. We will take the fol- 
lowing form 28 for the band structure in a cuprate supercon- 
ductor: 

^ q {eV) — to + f](cos q x a + cos q y a)/2 + t 2 cos q x cos q y 
+ £3 (cos 2q x a + cos 2q y a)/2 + t 5 cos 2q x a cos 2q y a 

+ £4 (cos 2q x a cos q y a + cos q x a cos 2q y a)/ '2 (5) 

where t = 0.1305, *i = -0.5951, t 2 = 0.1636, t 3 = 
-0.0519, U = -0.1117, t 5 = 0.051 (all in units of electron 
volts) and a ~ 4A. We will take a superconducting gap of the 
form A q = Ao(cosq x a — cosq y a)/2, and we will take the 
representative value Ao = 30 meV. 

For iron pnictide superconductors, we consider the 
BaFe2As2 family. We focus on the inner hole pocket located 
at kp ~ 0.1(7r/a), where again a « 4 A, with s-wave super- 
conducting gap A q = 12meV as in Bao.6Ko.4Fe2As2p^. We 



take the following simplified form for the energy dispersion 
of the inner hole pocket: 



£ q (eV) = 0.5(cos q x a + cos q y a) — 0.95 



(6) 



B. Tunneling between the Nanowire and the Superconductor 



Using the assumed form ip(x, y) — y ^ sm (i5r) f° r 

the electron wavefunction in the i th sub-band in the nanowire, 
we can write U(k x , q) in terms of the position space tunneling 
amplitude t(r, r') as follows: 



U{k x ,q) = 



dx J dy I d 2 r't(r,r')e lk ^ x -^ x 



, e-^' sin(^) 



W 



(7) 



Confinement of the wire in the y-direction destroys strict con- 
servation of y-momentum. As a result, we often encounter the 
following expression, which quantifies this non-conservation: 
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sm((q y w - ttj)/2) _ 
q y - irj/w 
? . sin((g y w + 7rj)/2) 
q y +irj/w 



(8) 



As w is increased, this function becomes more sharply peaked 
around q y = ±irj/w. In nanowires of width w w 100 nm, 
momentum non-conservation is small compared to the scale 
of the Fermi momentum in the superconductor since k F ~ 
ir/a m 250tt/w. 

We make further progress by considering several forms 
of the tunneling matrix elements t(r, r') at a clean inter- 
face between a higher- T c superconductor and a semiconductor 
nanowire, which we will refer to later: 

1. A smooth, uniform interface (with either a cuprate 
or pnictide) where t(r, r') = ~#( 2 )(r — r'). Then 
ti(k x ,q) = t\{k x ,q) with 



(9) 



2. A smooth uniform interface with a cuprate supercon- 
ductor such that electrons from the nanowire tunnel 
only into Cu 4s orbitals. Such a situation is pos- 
sible because these are the orbitals that extend fur- 
thest in the z-directions. Indeed, this is the dominant 
path for tunneling between copper-oxide planes in a bi- 
layer and also through insulating spacer layers between 
copper-oxide planes 30 . However, depending on the na- 
ture of the topmost layers of the cuprate superconduc- 
tor, electrons may tunnel instead into other orbitals as 
well. Therefore, this form of tunneling presupposes that 
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the topmost (presumably insulating) layers have been 
engineered (perhaps through molecular beam epitax- 
ial growth) so that electrons from the nanowire tunnel 
through the topmost layers and predominantly into Cu 
4s orbitals in the copper-oxide plane. The Cu 4s or- 
bitals hybridize with neighboring Cu 3d x 2_ y 2 orbitals 
to form a (nearly) half-filled band, so that 

t(r, r') = [6 W (r - r' - at) + <5 (2) (r - r' + at) 

Z7T 

- 5^ (r - r' - ay) - 6™ (r - r' + ay)] . 
ThenU(k x ,q) =t?(fc x ,q) with 
ti(k x ,q) = 2tS(k x - q x )g l (q y )(cos(q x a) -cos(q y a)). (10) 

3. A dirty or rough interface. When the interface between 
the nanowire and the superconductor is dirty or rough, 
momentum is not conserved during the tunneling pro- 
cess. For illustrative purposes, we consider the extreme 
case of i(k, q) = A(q), independent of k. The mis- 
match between the Fermi momenta of the nanowire and 
the superconductor no longer matters, but electrons in 
the nanowire are now coupled to different parts of the 
Brillouin zone, where the gap can have different signs. 

4. A nanowire on top of a step edge of a cuprate, as shown 
in Figure [T] We assume that the terraces are evenly 
spaced, with the terrace edges at x n — nl. We take 
I « 7nm, which corresponds to an angle 9 w 10°. In 
the clean limit, the tunneling amplitude is dominated by 
the terrace edges x n . We assume the steps are wide so 
that tunneling only occurs at x = x n when x' < x; as 
usual, there is no such restriction for y and y' . Then 
i(r,r') = ^J2„ dS ( x ~ x n)S(x cos9 - x')S{y - y'), 
where d is the length scale for the region in which tun- 
neling happens for each step. This yields ti(k x ,q) = 
i-(fc x ,q) with 

t s i( k x,q) = td^2 5 (k x - q x cos 6 + jQ)g l (q y ) (11) 
where Q = 2ir/l and j is summed over the integers. 




Cuprate Superconductor 



FIG. 1: Sketch of a interface between a semiconductor nanowire and 
a step-edge of a cuprate. 



C. Induced Superconductivity 

Integrating out the superconductor's degrees of freedom 
generates an effective action for the nanowire S c s = Snw + 
S' where 
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(13) 
(14) 
(15) 



It is important to pause at this point and consider these equa- 
tions. They are the same for s-wave and d-wave supercon- 
ductors, but there is a crucial difference in the latter case: the 
gap A q may vanish along certain directions in the Brillouin 
zone and, consequently, the single-particle energy in the su- 
perconductor y '^2 -f | A q \ 2 may vanish at certain points. At 
the nodal points on the Fermi surface of a d x i_ y i supercon- 
ductor, the denominators in Eqs. T3fT5 vanish quadratically 
at uj m = 0. Therefore, these integrals will diverge loga- 
rithmically unless the numerators also vanish. However, for 
an infinitely-long nanowire, momentum conservation along 
the wire (or momentum conservation up to a multiple of Q) 
prevents any coupling between low-energy electrons and the 
nodal points of the superconductor. Consequently, the numer- 
ators in Eqs. T3p5 are zero at the nodal points, where the 
denominators are dangerous. However, in finite-length wires, 
there will be a coupling to the nodal points, and we study this 
coupling perturbatively in Section[V] 

From the total action S e g — Snw + S', we obtain the 
spectrum from the poles of the Green function: 



cr 



( k r 



2m* w 2 



Sijdk^k'^ 



ak x a y r z + A' k Kij (u)T x (16) 



Here, we have analytically continued iu: n 



iS and writ- 



ten, e.g. r fcx fe' 



r k x k' 



id). The smallest positive 



pole of this equation is the gap. 

If momentum is conserved in the x-direction, then r, e', A' 
are all diagonal in k x . If we further assume the eigenvalues uj 
are much smaller than A q , then we can drop the ui dependence 
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of nfe^jH, and A^^). Finding the poles of 

Eq. ( fl6| reduces to finding the eigenvalues of the matrix M: 



M = (8ij + rij(k x )) 1 (ak x a y T z + A l ij (k x )T x 



( 



fi 2 (7Ti) 2 

2m*w 2 



Since w 3> a, momentum non-conservation in the y-direction 
is small on the scale of the Fermi momentum of the supercon- 
ductor. We make the following approximation, for the mo- 
ment, in order to take a qualitative look at the induced gap 
(we will use the full expressions when we turn to a more care- 
ful computation of the induced gap): 



9j(ly) 



Aj+l p -iq y w/2 



2i\/w 



-[S(q y -wj /w)-(-l) 3 S(q y +TTj /w)]. 

(18) 



In this limit, we can simplify the expressions in the previous 
subsection for t u and t n . Now, the effective action S s s takes 
the form: 



<S e ff = 



2tt 

J 



mA 



' l e fcx r 2m*w 2 ~~ V ~ 6 k x A,m) T z 

+ V x a x - A^,^^]^^^ 



with Eqs. ( 13 1-( 15 1 now taking the simpler form 
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(20) 
(21) 
(22) 



We now see that, in the limit in which we replace gj (q y ) by a 
sum of (5-functions, there is no coupling between the nanowire 
and the nodal points in the superconductor for generic val- 
ues of k x (including the expected Fermi momenta for realistic 
nanowires). In these equations, t{ kx is given by 
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ti(k x ,q)ti{k' x ,qy 



1*1 fort,; 



tf 



~ \t\ 2 (cos k x a — cos k v a) 2 for ti = t\ 
The induced superconducting gap function is: 
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At the Fermi surface, in the static limit, this is 
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\U,k F \ +£kF,7ri/w \^k F ,7ri/w\ 



(26) 

Note that £ kp ( may not vanish due to the mismatch between 
the Fermi momentum of the nanowire and that of the super- 
conductor. This mismatch is one of the limiting factors for 
induced superconductivity. 

From the single-particle spectrum obtained from Eq. (17) , 
we see that the single-particle gap at kp in the i th sub-band 
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(27) 



where Eso = okp^ x and kF,x is the Fermi momentum in the 
x-direction. 



D. Renormalized Parameters and Topological Phase 
Transition 



The effective action S e s takes the form of the action of 
a nanowire with an applied pair field A' k im , as may be 
seen in the w a limit in Eq. (19i. However, the pa- 



rameters fi, a, e, V x are renormalized compared to those of 
an isolated nanowire as a result of the coupling to the su- 
perconductor. In the w 3> a limit, this is simply renormal- 
ization by a factor of (1 + r) _1 , as may be seen by making 
the change of variables ^ m ,i,K -> ^m^fc, (1 + 



^m,i,k x -> $m,i,k x 0- + r k x ,i,m) 1/2 in E q- ( |19) , which now 
takes the form: 



2^. ^ r 



+ ( e fe! d - Mfe^i.m) T z + (Vx) 1 k^i, m (7 x-^k n ^i, m T x]^rn,i,k :c 



(28) 



where 



{V x f^i,m = V x /(1 + r kxA , m ), al\ m = a/(l + r k ^ m ), 

/'/':'',.„ = (M + rk.,i, m ). (29) 

In particular, the Zeeman splitting V x takes a renormalized 
value. Since V x is independent of \t\ 2 but r is proportional 
to \t\ 2 , larger tunneling effectively reduces the g-factor of the 
nanowire. As we will see in subsequent sections, the condi- 
tion to be in a topological phase sets a lower bound on V x nd . 
Hence, a reduced effective g-factor would require a larger 
magnetic field to reach the topological regime. One possible 
concern is that the larger required magnetic field would de- 
stroy superconductivity, but in the high-T c materials that we 
consider, the critical field is large enough that this is not the 
case. Even if we could induce a superconducting gap in the 
nanowire on the order of the gap in a cuprate superconductor, 
A md ~30 meV, the threshold to be in the topological phase 
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would require V x nd to be at least this large, corresponding to 
an applied field of about 40T for the values of r ~ 1 that we 
consider in this paper. While this field is very large, it is still 
less than H C 2 in a cuprate superconductor. Moreover, as we 
will see, the induced superconducting gap is typically a few 
meV and, therefore, easily achievable with magnetic fields of 
a few Tesla, which will have negligible effect on a high-T c 
superconductor. 

When searching for the topological superconducting phase, 
it is the renormalized parameters and not their counterparts 
measured before the wire is coupled to the superconductor 
that determine the onset of the phase. The topological phase 
is characterized by the existence of Majorana zero modes at 
the end of the wire, which occurs when a condition of the fol- 
lowing form is satisfied: 



where, 



(30) 



The parameters that enter this condition are renormalized. 



E. Majorana Zero Modes 

We now consider the possibility of having Majorana zero 
modes at the end of the wire. For simplicity, we will assume 
in this subsection that the Fermi level lies in the lowest band 
and take sub-band index i = 1. We will suppress the sub-band 
index and write ip = ipi = i. The Hamiltonian of the nanowire 
can be written in the following form in real space, where all 
parameters now correspond to their induced values after cou- 
pling to the superconductor, described in the previous section 
(we drop the superscripts 'ind' in order to avoid clutter): 



dx 



H 



We assume that the nanowire lies along the negative cc-axis 
and terminates at x = 0. This condition can be realized by 
setting [i(x < 0) = fio, and fi(x > 0) = — oo. 

Note that such a BdG Hamiltonian is justified in the d x i_ y i - 
wave case if the nanowire is decoupled from the nodal points 
of the superconductor. This is the case in the limit of an 
infinitely-long wire, as we saw in Section |ll"C| However, when 
we consider Majorana zero modes at the ends of (the topolog- 
ical portion of) a wire, we are necessarily faced with a sit- 
uation in which momentum along the wire is not conserved, 
so there will be some coupling between the zero modes and 
nodal excitations. We set this coupling to zero here and treat 
it perturbatively in SectionfV] 

Then our Hamiltonian in the Nambu basis ^(x) = 
(ip].(x), ipj_(x), ipi(x),i/j^-(x)) can be written as, 



H = J dx^{x)H BdG ^{x) 



H 



BdG 



dx 



8(x - x') (^- t y^f - n{x) + V x a x ^j 



+ 5(x — x')iaa y d x + A(cc, x')a z T x 
which gives the following BdG equation for E = 0: 



(32) 



BdG ' 



{ut( x )' u i( x )> v i( x )' v t( x )) T = (33) 



Since the BdG Hamiltonian is real, we can have real solutions 
for Majorana zero modes. After imposing particle-hole sym- 
metry for a real solution, we can set /i(x) = Xu^ n {x) with 
A = ±1. The BdG equation for E = can be written as, 



dx' 



n-di 



-s(x-x')i-^f+n 
V-(x,x') 



V+(x,x') 



u t (x') 



(34) 



where 



V± (x, x') = V x S(x -x')±X A(x, x') ± aS(x - x')d x (35) 

with 3 constraints: [u^^(x — 0)] = and normalization. 
Assuming u^n(x < 0) oc e zx , the existence of a zero mode 
requires at least three roots Zi with positive real part, so that it 
is normalized and localized at the end x = 0. 



{^U x ) \-%^ - K x ) + ioca y d x + V x a x ^j i/> a >(x) 
dxdx'[{A{x,x')i)\{x)^\{x') + h.c.]\, (31) 



III. SUPERCONDUCTORS WITH s± PAIRING 
SYMMETRY 

Recent results on iron pnictide superconductors suggest 
that its superconducting order parameter has s± pairing sym- 
metry. The Fermi surface has several components, and the 
sign and size of the gap vary from one Fermi surface compo- 
nent to another but the gap does not change sign as any Fermi 
surface component is encircled. The smaller hole pocket cen- 
tered at the r point has a gap A w 12meV (see, e.g. Ref. 
l29t . Therefore, in the clean interface limit with momentum 
conserving tunneling, we can induce a large superconduct- 
ing gap at the Fermi points of a nanowire by bringing the 
Fermi momentum of a nanowire near to this hole pocket's 
Fermi surface. The Fermi level for the smaller hole pocket 
near the T point ranges from 0.1 to 0.3 {it /a) depending on 
the doping level 32 . For a = 4A and either w = 50 or lOOnm 
with 5 sub-bands occupied, there is a mismatch between the 
Fermi momentum in the nanowire, which is about 0.02 to 0.04 
(n/a), and the Fermi momentum in the pnictide. This reduces 
the induced paring potential, as we will see in detail below. 
For simplicity, we neglect the momentum dependence of the 
gap in the inner hole pocket, Afc = Ao = 12meV, as ob- 
served in Ref. |29]and assume momentum-conserving tunnel- 
ing, t u (k,q), as introduced in Section II B We see from Eq. 
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FIG. 2: Fermi surface of a typical iron based superconductor. Two 
hole pockets at the center (kp = 0.1 ~ 0.3(-7r/a)) and particle pock- 
ets at the corners have superconducting order parameters of opposite 
sign. There are five sub-bands in a nanowire with w — 50nm and 
kp = 0.04. The 10 lines inside the dashed circle correspond to 5 
sub-bands in the nanowire and their endpoints are the Fermi points 
of these sub-bands. 
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(36) 
(37) 



We find that Tk x ,i,j and A' k i ■ are essentially momentum 
independent and diagonal in sub-band indices i, j (see Ap- 
pendix A), so we can assume that the induced gap behaves 
like an s-wave superconducting gap. Using Eq. (34i with 
A(x,x') = Af Ld and /irj = Mi> we 8 et a quartic equation for 
z for each sub-band: 



1 



z 4 + (/Zj + a 2 ) 



~2\ 2 



2AA,5z + /I, 2 - Vl + A 2 = 0, (38) 



where x = 



m (V.f; 

h 2 ■ 



' '■ m ot 



F, = and A. t = 

In this case, the condition for 



r, 2 



A\ nd with u tj ; cx e 2 



the i th sub-band to be in the topological phase is given by-^ 
/I; 2 — V x + Af < 0. Of course, our main concern is that the 
highest occupied sub-band (or an odd number of sub-bands) 
be in the topological phase. 

We now compute the induced pairing potentials for the pa- 
rameters described above. We consider aw — 50nm wire 
where i — 1,2,3,4 sub-bands are occupied, and only one 
chirality-split branch of the i = 5 sub-band is occupied. We 
set /xg = 5meV and V x — 15meV to ensure the topologi- 
cal phase only in the 5th sub-band. For t = 60meV, we get 
r*a ~ 0.12 which does not renormalize the effective action 
much. Then the diagonal element of induced gap in each sub- 



band is (in meV): 



a ind 
^1 


— 1 ^7 


a ind 
^2 


= 1.37 


a ind 
^3 


= 1.37 


A!t nd 


= 1.38 


A i„d 


= 1.38 



(39) 
(40) 
(41) 
(42) 
(43) 



The topological phase transition occurs at V x = 
\J /i 2 + | A" ld | 2 . Therefore, the magnetic field required to ob- 
serve the topological phase transition is a few tesla if g ~ 50. 
Since the iron pnictide superconductors have large H C 2 ~ 
50T, the induced gap will not be significantly suppressed by 
the applied magnetic field in the region of phase transition. 

Although the induced pairing is large, as may be seen from 
Eq. ( 39 1, the resulting single-particle gap, which is the most 



physically-relevant quantity, is significantly smaller according 
to Eq. (|27). For V x = 15 meV, A^ = 0.16 meV. 

The other types of interfaces are less interesting for pnic- 
tide superconductors. At a dirty interface, electrons from the 
nanowire will be able to tunnel to all of the components of the 
Fermi surface, which will suppress the induced gap since the 
electron pockets are expected to have superconducting gaps of 
the opposite sign. A step-edge interface could, similarly, al- 
low tunneling to the electron pocket, which will suppress the 
induced gap. 



IV. SUPERCONDUCTORS WITH d x 
SYMMETRY 



,2 PAIRING 



A. Clean, Uniform Interface between a Cuprate 
Superconductor and a Semiconductor Nanowire 

In this section we consider superconductors with d x 2_ y 2 
pairing symmetry, such as the hole-doped cuprates. Then 
A q = ^ (cosq x a — cos q y a). First, let us assume that the 
interface between the SC and the NW is clean and flat and, 
therefore, momentum conserving. The direction of the NW is 
assumed to be parallel to the a-axis of the SC. For the simplest 
possible form, ti(k x ,q) = t^(k x ,q), we see that 



dq y 



1*1 9i(qy)9*i(q y 



(2tt)2 ( 



•a. 



Ai 



(44) 



and 



a; 



l*l 2 9i(q y )9j(q y ) x 



(cos k x a — cos q y a) 



A 



(45) 



kx,q y | 



Since 9i(q y ) is peaked at ±iri/w, the momenta {k x , q y ) will 
be far from the Fermi surface of the superconductor for the 
first few nanowire sub-bands. Therefore, A',. - ■ „ will be 
strongly suppressed. Focussing, again, on the case of four 
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and a half occupied sub-bands with /15 = lmeV and V x = 
2meV in the nanowire with w = 50nm, we find r^^j has 
very weak dependence on k x and is almost diagonal in sub- 
band index. We obtain r, ~ 0.47 for t = leV. We also find 



0.05 



> 



-0.05- 




k (xl0°m 



FIG. 3: Induced gap as a function of k x in each sub-band. Dots 
denote the integration values. Plots of the k x dependence of a d- 
wave gap (with k y fixed to the value expected for the corresponding 
sub-band) are presented for guidance. 

that the induced gap for each sub-band has the k x dependence 
expected of a d-wave SC with k y fixed to the value expected 
for each sub-band, as one can see in Fig. [3] The value of the 
induced gap at the Fermi momentum of each sub-band is (in 
meV) 



a ind 


= -0.052 


(46) 


a ind 


= -0.038 


(47) 


kf,3 


= -0.014 


(48) 


a ind 


= 0.018 


(49) 


Ai nd . 


= 0.059. 


(50) 



For V x =2 meV, this gives A qp = 0.02. 

We now consider Majorana zero modes at the end of 
such a nanowire. If we make the approximation g(q y ) ~ 
5 (q y - ™) + 6 (q y + ^) following Eq. we can take 

A(x,x') w W (a 2 d 2 x - (^) 2 ). The BdG equation can 
then be written as in Eq. ([34| with: 



V± = V X ± XWi 



[a 



2 d: 



(51) 



The above BdG equation can be written as a quartic equa- 
tion for z with real coefficients. 



A' 



z 4 + 2AAz 
2tt 2 AA 



1 + + 



2tt 2 A 2 



z + ft 2 



V 2 



A 2 7T 4 



0, (52) 



where x = ro fc " z , /i 



ft fJ-o 



V x = 



ft V x 




FIG. 4: Fermi surface of a cuprate superconductor. The 10 lines 
inside the dashed circle correspond to 5 sub-bands in the nanowire. 
The ends of the lines are the Fermi points of the nanowire, which are 
located far from the Fermi surface of the SC, as may be seen from 
the figure. 



channel, — z^s are the solutions for A = — 1 channel. Since 
the coefficients are real, if zi is solution, z* is also a solution 
for same channel. 

1 . When ft 2 — V 2 + A ,J < 0, there is at least one negative 
real root and one positive real root. Also, the product 



of the four roots Z; is 



< 0. When all 



1+4A 2 

roots are real, we have three positive roots for either 
A = 1 or A = —1. When two of the roots are complex, 
then the four roots can be written z\ > 0, z 2 < 0, 
23 = a + bi and Z4 = a — hi, and we again have three 
roots with positive real part for either A = 1 or A = — 1. 
Therefore, we have a MZM in this case. 



2. When fi 2 - V 3 



> and all four roots are real, 



and w — ^jrr-w. Note that when z^s are the roots for A = 1 



there are two different cases. When two of them are 
positive and two of them are negative, we do not have a 
localized solution for zero energy. When all four roots 
have same sign (which is positive for either A = 1 or 
A = — 1), we have two MZMs at the end of nanowire. 
However these two localized states are at the same end, 
and they will split into two states with E > and E < 
by interaction. 

3. When p? — V 2 + ^-Z > 0, and two roots are complex, 
the other two roots, if they are real, will have same sign 
since Yit=i z i > 0- ^ tr, e other two roots are also com- 
plex, those two also have the same real part. Then it is 
similar to case 2. 

We do not consider the situation in which the above equa- 
tion has a double root or two purely imaginary solutions in 
case 1 since those cases are sets of measure zero in parameter 
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space. Therefore, we can conclude that there is a single Majo- 
rana bound state with zero energy when jl 2 - V 2 + < 0. 
Indeed, we can extend this analysis to the nth sub-band where 
k y = nir/w, and the condition for topological phase is given 

„„ r.2 T/2 , A 2 K 4 7r 4 s r, 



B. Tunneling Through Cu 4s Orbitals 

We also consider the case in which electrons tunnel from 
the nanowire only into Cu 4s orbitals. These orbitals hy- 
bridize with neighboring Cu 3d x 2_ y a orbitals to form a half- 
filled band at the Fermi energy. Then ti(k x ,q) — t1(k x ,q). 
We now have the intriguing form 



7^2 \t? 9i{%)9*j{% 



^7, 



(cos k x a — cos q y ay 



(53) 



Although this expression is more sharply peaked at the anti- 
nodes, where the gap is large, it suffers from the same limita- 
tion as the previous one, namely that it is strongly suppressed 
by the Fermi momentum mismatch between the nanowire 
and the superconductor. This is exacerbated by the fact that 
nanowire electrons have momenta that are close to k = (0, 0) 
(on the scale of the cuprate's Fermi momentum), and both 
i(k, q) and A^ are very small there. So, although this type of 
tunneling manages to avoid the nodes, it does not really allow 
the nanowire electrons to tunnel to the superconductor's anti- 
nodes. Consequently, the induced gaps are very small, as we 



show in more detail in Appendix A 3 



C. Dirty or Rough Interface 



We now consider the case of a dirty or rough interface. In 



the extreme case introduced in Section|HB] Eq. ( 15 i becomes 



d\ |A(q)| 2 A q 



(2-) 2 aA 



|A„ 



The right-hand-side is independent of k. Moreover, since A q 
is odd under rotation by it/ 2 while the rest of the integrand is 
even, the right-hand-side vanishes after integration, and there 
will be no induced gap. 



D. Nanowire - Step Edge Interface 

We now analyze the situation in which the nanowire is on 
top of a step edge surface of a cuprate, as shown in Figure [T] 
For simplicity, we assume that the terraces are evenly spaced 
so that the terrace edges are at x n — nl. In the clean limit, 
the tunneling amplitude is dominant at the terrace edges, x n . 
We assume the tunneling matrix element ij (k x , q) = t\{k X: q) 



discussed in Sec II. For small angle 9 ~ 10°, we can approxi- 
mate cos 9 ~ 1, and 



A' 



E 



d 2 q 2 

j— -2 \t\ 2 8(k x -q x + mQ) 



x5(k x -q x + n 2 Q)— — —2 (54) 

^m + £g + |Ag| 

At the densities under consideration in the NW, kp << Q, 
and we can ignore the contributions from k x = k x + (ri2 — 
n i)Q w i m n i n i- By integrating over k' x as we did before, 
we get 



a; 



E 



Afc x +raQ,g y 

w m + ^k x +nQ,q y + \^-k x +nQ,q y \ 



(55) 



and 



r k,. 



E 



dq y 



\t\ 9i(Q y )9j{q y ) 



+ \ Ak T +nQ,q v | 



(56) 



with k x + nQ G first B.Z. 

Again, in the limit of small inter-band coupling, the induced 
gap for each sub-band is given by 



A ITid 



1 + nc F ,i,o 



(57) 



It is now possible for (k x +nQ, q y ) to lie on the Fermi surface 
of the cuprate superconductor, so the induced gap is no longer 
suppressed. 

To estimate the size of the induced gap in each band, we 
take /is = lOmeV and V x = 30meV for a nanowire of width 
w = 50nm. Again we find r^^.j is constant in k x and diag- 
onal in (i,j). Choosing t = 25meV gives rj ~ 0.11 to 0.13. 
We find the induced gap for each sub-band has a very dif- 
ferent k x dependence than in the uniform tunneling case: the 
induced gap is only weakly-dependent on k x and, in particu- 
lar, it does not change sign. The induced gap in each sub-band 
at its Fermi momentum is 



a ind 
^k F ,l 

a ind 
^■k F ,2 

a ind 
A ind 



A 



ind 



-3.04 
-3.09 
-3.15 
-3.20 
-2.64, 



(58) 
(59) 
(60) 
(61) 
(62) 



in meV. For V x = 30 meV, this gives A qp = 0.2. In reality, 
tunneling will not be perfectly momentum-conserving modulo 
Q. However, the basic result should still be valid: if momen- 
tum non-conservation is much larger in one direction than the 
other then a large gap can be induced. 
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FIG. 5: Fermi surface of an interface between a semiconductor 
nanowire and a step-edge of a cuprate. Now induced superconduct- 
ing gap gets contributions from the colored area. 




-2 



-1 



-a-i=l 
-*-i=2 
H-i=3 

-0-i=5 
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FIG. 6: Induced gap as a function of k x in each sub-band. Dots 
denote the integration values. The lines are piecewise linear interpo- 
lations of the points. The step occurs when one extra multiple of Q 
contributes. Note the contrast with Fig. [3] 



We now consider Majorana zero modes at the end of the 



wire. As discussed in Appendix A 4 we find the gap to vary 
very slowly with k x since the biggest contribution comes from 
k x + jQ ~ tt/o where A/, and ^ does not change rapidly 
with k x . Then, we are effectively in the case of an s-wave 
superconductor as we were in Sec III, and the characteristic 
equation for z takes the form: 

* z 4 + (ft + a 2 ) z 2 - 2AA,<5z + ft 2 - V 2 + A 2 = 0. (63) 

By essentially the same analysis as above, there are three so- 
lutions with positive real part so long as jl 2 — V 2 + Af < 0. 



V. MIXING BETWEEN THE MAJORANA ZERO MODE 
AND GAPLESS BULK EXCITATIONS IN A d x 2_ y 2 -WAVE SC 

High-Tc superconductors with d-wave pairing symmetry 
have four gapless nodes in the two-dimensional Brillouin 
zone. Thus far, we have ignored these low-energy excitations 
because the tunneling matrix elements t(k x ,q) that we used 
did not couple electrons in the nanowire to these nodal excita- 
tions so long as the wire was infinitely-long. In a finite-length 
wire, however, momentum along the wire is not conserved, 
and there will be some coupling between low-energy electrons 
in the wire - especially the zero modes - and nodal excitations 
in the superconductor. In addition, impurities could scatter 
nanowire electrons to the nodal points. In this section we add a 
tunneling term to the action coupling the Majorana zero mode 
to the superconductor. We then calculate the self-energy of 
the Majorana mode perturbatively to analyze whether it will 
survive or decay into the bulk. 

The tunneling action coupling the j th sub-band Majorana 
zero mode localized at the end of the wire at x — to the 
fermionic excitations in the superconductor is 



/0 rw p 

dx dy dVu(r,r')7m,j(i-) 
,„ -oo JO J 



= / %,-(k)7j(wm) [cZ. 



m , — k 



-m,— kj 



(64) 



where v(r, r') is the tunneling amplitude from the Majorana 
zero mode at r to a state at r' in the superconductor. We take 
v(r,r') to be real. On the superconductor operators we have 
explicitly written the spin superscript u. In going from the first 
to the second equality, we have assumed a simplified form for 
the real Majorana zero mode: 



,r) = 2 



( T^y \ 



sin I — I Q{-x)e zx l3 {uj m ) (65) 
If \ w J 



We use the Fourier convention for v lj (k) 

v l3 (k) = 2^-J dx J dy J dV w(ry) (66) 

( ^3V\ c zx+tk-r' 



x sm [ — — I e~ 
\ w I 



(67) 



Then the self energy S 7j . 7j .(a; m ) = (jj(oJ m )jj{-u} m )) is 
given by: 



S 7 ~ (w m ) = / d 2 k |u 7i (k) 



, d 2 k l u 7 3 ( k )| {A k +iuj m ) 



(68) 



For simplicity, we take to be of the form ^ = 1 1 (cos k x a - 
cosfcj,a)/2 + ti cos k x a cos k y a and Ak = Aq(cos k x a - 
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cos k y a). We consider three types of tunneling correspond- 
ing to the three interfaces described in Sec. |IIB| for which the 
tunneling matrix elements are given explicitly by 



can be handled together: 



ImS^(w) 



<( k ) = 
<-(k) = 



t\/zgj(-ky) 

z ~\~ 

h£MzM 2 (cosk x a- cos k y a) (70) 
tdyfzgj(-k y ) 

~y g— (z-\-ik x cos0)l 



V% Aeko \J c\p 2 + c\q 2 



(78) 



(69) 



(71) 



A simple change of variables yields the decay rate T(uS) to 
leading order: 



(79) 



A. Zero Temperature Self Energy 

We obtain the retarded self-energy by taking iu) m —> uj + irj 
and then use the identity lim n ^.o+ ( x +i V ) = ^(x) ~ ^ % ^{ x ) 
to obtain 



where a and j3 are dimensionless numbers, and s — a 2 z or 
d 2 z in the cases of uniform and step-edge tunneling, respec- 
tively. In the case of tunneling through the Cu 4s orbital, we 
have: 



v2 ^kek 



y/cfp 2 + c$q- 



4 «2 



ImE^H = 2 / K(k)|> + A fc )|V fe (4 2 + A£)| 

Aek 



so that 



where k satisfies u 2 — £^ + A 2 ^, and 
ReS7 7 (o;) = lira 2 / d 2 k |w 7j (k 



(72) 



T n (uj) oc 



a 2 z\ u J i 



(80) 



(81) 



2 (A fc + c)(4 2 + A 2 fc -^) 
(73) 

From here on, Afc disappears from the numerator because it 
is odd under exchange of k x and k y , while all other terms 
are even. We now explicitly calculate the real and imaginary 
parts. 



1. Imaginary Part of Self Energy 



In all cases, the decay rate is suppressed by a small coeffi- 
cient. The decay rate for tunneling through the Cu 4s orbital 
is further suppressed at low energies by its cubic dependence 



2. Real Part of Self Energy 

The real part of the self-energy follows similarly: 

w(c\p 2 + c 2 q 2 - uj 2 ) 



ReYP r ~<{u]) = lim 2 / dpdq\v u ' 



J ' ' {c\p 2 + C 2 q 2 — UJ 2 ) 2 + T] 2 



For small ui, the dominant contribution to Eq. ( 72 1 comes 
from momenta near the nodes, which we denote by overbars: 
(ifoj, iky). We expand the momenta around the nodal point 
(kx, ky) as (k x , k y ) = (k x + p + q,k y + p — q) and expand 
similarly around the three other nodal points. Expanding 
and Afe about the nodal points yields ^ = c\p and Afc = c%q, 
where c\ = —t\asm.k x a — 2i; 2 a sin k x a cos k y a and c 2 = 
2Ao<x sin fc x a. Now the condition k e ko is given by u 2 — 
c\p 2 + c 2 q 2 . We also linearize the tunneling strengths: 



ClC 2 



(82) 



<.(k) = hlM-M =v u 



z ~\~ xkx 



(74) 



<,( k ) 



tyfetoi n y\ qasili k xa = v n qa (75) 



<(k) 



z ~\~ %k x 
td\fzgA—ky 



where A is a high energy cut-off. Hence, in the small u limit, 
the real part of the self-energy is logarithmically-divergent. 
Therefore, the would-be pole at zero energy in the zero-mode 
Green function has zero weight. This is a sign that the zero 
mode does not survive the coupling to nodal quasiparticles. 
Eventually, the zero mode will leak into the bulk of the d x 2_ y 2 
superconductor. However, the divergence is only logarithmic 
because there is little phase space at the nodes, so this leakage 
occurs very slowly. Moreover, the coefficient in front of this 
divergence is small for the reasons noted above. 
For tunneling through the Cu 4s orbital: 



7j v 1 -y g — (z-\-ik x cos 9)1 

At two of the nodal points, q is replaced by p in u™. . Note that 

\v u ' n \ 2 (xa 4 z/w, \v s \ 2 <xa 2 d 2 z/w (77) 

are small because 100 < w/a < 250. In addition, d ~ a and 
1/z ^> a, which further suppresses these tunneling parame- 
ters. The cases of uniform tunneling and a step-edge interface 



(76) JteSZPYw) = lim / dpdq- 
n^-o J 



2\v n \ 2 a 2 uj 



\ci4\ 



2\v n \'a 2 uJC 2 q 2 {c 2 p 2 + c 2 q 2 - uj 2 ) 
[c\p 2 + c\q 2 — uj 2 ) 2 + i] 2 

uj 2 ln^/VA 2 -uj 2 )] (83) 



\ c ' 
A 2 
2 



Hence, in the small uj limit, the real part of the self-energy 
goes like w. There is no divergence because tunneling gets 
weaker as the nodes are approached. Therefore, there is a 
pole at zero energy and the zero mode survives. 
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B. Lifetime at Non-Zero Temperature 

At non-zero temperature, the preceding calculation is mod- 
ified in two ways. First, the Matsubara frequencies must take 
values £j m = (2m + 1)tt/@. However, so long as we ne- 
glect electron-electron interactions in the nanowire, there is 
no sum over Matsubara frequencies; the Matsubara frequency 
of the zero mode is precisely the same as that of electrons 
or holes in the superconductor. Therefore, this manifestation 
of non-zero temperature makes no difference. The second 
way in which non-zero temperature can enter is through the 
Green function of nodal fermions in the high-T c superconduc- 
tor. At non-zero temperature, these fermions have a lifetime 
proportional to T (see, for instance Ref. [33l . Therefore, the 
superconductor's Green functions are modified according to 
lj -> lj - Ssc(w,T) with ImY, sc (u,T) w T for w < T. 
Consequently, we now find that, for lj < T: 

T U ' S {T) ~ TlnT 
r"(T) - T 

The coefficients depend on the tunneling strengths v in the 
same way as at zero temperature, so the decay rates will again 
be suppressed by numerical factors. 



smaller systems or higher temperatures. However, the spin- 
orbit coupling strength is equally, if not more, important for 
determining the energy and length scales that protects the zero 
modes. We have seen that A qp is consistently an order of 
magnitude smaller than A lnd ; this occurs because the spin- 
orbit energy is an order of magnitude smaller than the Zeeman 
energy which, in turn, must be large because the induced su- 
perconducting gap is large. Therefore, it would be interesting 
to see how the choice of interface geometry - for instance, the 
terrace structure that helps a d x 2_ y 2 -superconductor induce 
superconductivity in a nanowire - can help enhance spin-orbit 
coupling in the nanowire. 
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Appendix A: Numerical Values of the Gap and Renormalized 
Dispersion 



VI. DISCUSSION 

Higher- T c superconductors are not a panacea for enhancing 
the stability of Majorana zero modes. Part of the beauty of 
nanowire-superconductor proposals for Majorana zero modes 
is that they are relatively insensitive to the details of the su- 
perconductor. However, since the cuprates and pnictides (and, 
perhaps, all higher- T c superconductors) have gap functions 
that change sign in the Brillouin zone, their details necessarily 
matter. Indeed one might, initially, expect that, as a result, it 
is impossible to use them to induce robust topological super- 
conductivity in a nanowire. However, we have seen that it is 
possible, with the right type of interface, to do so. In fact, for 
a nanowire in contact with a cuprate superconductor at a step- 
edge surface, a pairing gap of ps 3 meV can be induced, and 
it is conceivable that a pairing gap as large as 15 meV could 
be induced with sufficiently strong tunneling at the interface. 
Even in the case of a pnictide superconductor, a gap of « 1 
meV seems achievable. 

In the case of the cuprates, there is a second difficulty: the 
presence of nodal fermions. Naively, these should immedi- 
ately wipe out the possibility of topological superconductivity 
and Majorana zero modes. However, the nanowire's electrons 
couple weakly to the nodal points, so the decay of the zero 
modes may be slow. They will always decay into the bulk 
(except in the special case in which tunneling occurs solely 
through the Cu 4s orbital), but for w /a sufficiently large, this 
may not be a quantitatively larger effect than the decay into 
impurity states in an s-wave superconductor. 

The primary focus in this work has been to increase the en- 
ergy and temperature scale associated with superconductivity. 
This could help make Majorana zero modes more robust in 



We present numerical values for the matrices r, e' and A' 
defined in eqns [T3p4 evaluated at the Fermi momenta of the 
highest filled band in several different cases. From now on, the 
units are (t/leV) 2 for r and t 2 /(leV) for e' and A'. Values 
for the induced gap A md are presented in plots as a function 
of t. 



1. Pnictides: s± pairing 

We assume that the Fermi level lies between the chirality 
split bands of the fifth subband and take the representative 
values /15 = 5meV and V x = 15meV. The Fermi momentum 
of the fifth band is kp^ = .086nm _1 . For a wire of width 
50nm, we find the matrices at ui m = to be: 



/30.6 .036 .062 \ 
31.2 .098 1 
036 32.3 .190 
.098 33.78 
V.062 .190 35.9/ 



/ 1.52 -.0003 -.0005 \ 

1.53 -.0008 ' 

.0003 1.55 -.0016 

-.0008 1.58 

V-.0005 -.0016 1.62 / 



Since the pairing function is constant the matrix A' is propor- 
tional to the matrix r. The matrices evaluated at the Fermi mo- 
menta of the other bands yield similar values. Hence we con- 
clude that r and e' are diagonal in both k x and sub-band index 
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(the latter follows from the small values of the off-diagonal 
matrix elements above). Consequently, we can find the in- 



duced gap for each band from Eq. ( 25 1 as a function of the 
tunneling strength t (defined in Eq. (|9j)), shown in Fig [7] The 
figure also shows the induced gap for a wider wire of width 
w = lOOnm, for which the matrices r and e' are similar to 
the above matrices, but with smaller off-diagonal elements 
because y-momentum is more nearly a conserved quantity. 
Since r, e', and A' remain small even for t as large as lOOmeV, 
it is possible that even such a large tunneling matrix element 
can be treated as a small perturbation. However, even for more 
a more modest tunneling matrix element such as 50meV, the 
induced gap is approximately 1 meV. 
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FIG. 7: Induced gap for a cuprate superconductor with step edge 
tunneling (top lines) and an s± superconductor (bottom two lines) as 
a function of tunneling strength. Within each pair the larger induced 
gap corresponds to the narrower wire. 
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Similar matrices are obtained at the Fermi momenta of the 
other filled bands, although there is variation in A' with k x 
which can be positive or negative, but is always very small. 
The r and e' matrices are effectively diagonal, although A' is 
not. However, since the elements of A' are so much smaller 
than the energy splitting between bands, we can still treat the 



effective action in Eq. ( 19 1 to be diagonal in bands to estimate 
the induced gap from Eq. d25}. The induced gap is extremely 
small; we plot it as a function of the tunneling strength t 
in Fig [8] along with the induced gap for a wider wire with 
w = lOOnm. 
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FIG. 8: Induced gap for a cuprate superconductor with uniform 
tunneling as a function of tunneling strength for wires of width 
w = 50nm (top) and w = lOOnm (bottom) 



2. Cuprates: uniform tunneling 



We assume again that the Fermi level lies between the chi- 
rality split bands of the fifth subband and take the representa- 
tive values /i 5 = lmeV, V x = 2meV. The Fermi momentum 
of the fifth band is kp^ = .034nm~ 1 . For a wire of width 
w = 50nm we find the following matrices at LJ m = 0: 



.465 


.00002 





.466 


.00004 



\ .00003 



,00002 


.467 


,00008 



.00003 \ 
.00004 ' 

.00008 
.469 
.470 



3. Cuprates: tunneling through Cu 4s orbitals 

Using the same representative values fi 5 = lmeV and 
V x = 2meV we can find the matrix elements for tunneling 
into Cu 4s orbitals described in Sec III B I These elements are 
greatly suppressed by the factor (cos(q x a) — cos(q y a)) 2 in 
Eq. (10 1 because of the momentum mismatch between the 



wire and the superconductor and are two or more orders of 
magnitude less than the uniform tunneling described in the 
previous section. Hence, we will not discuss this case in any 
more detail. 
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4. Cuprates: tunneling through the step edge interface 



Following Sec IV D we anticipate a sizable induced gap in 
this case because k x = q x only modulo Q, so we choose a 
larger value V x — 30meV to accommodate this. We choose 
/i5 = lOmeV. We again assume that the chemical potential 
lies between the chirality split levels of the fifth subband, 
which has Fermi momentum kp 5 = .12nm . We calculate 
matrix elements evaluated at this momentum and at u m = 0: 



/ 169 .050 .085 \ 
' 169 .126 ' 



050 172 .254 
.126 176 
V.085 .254 180 / 



/-10.8 .001 .002 \ 
-10.8 .003 ] 
.001 -10.9 .005 
.003 -11.0 
V .001 .005 -11.1 / 



At 



-4.37 


-.0005 





-4.41 


-.001 



\ -.0008 



-.0005 


-4.46 


-.002 



-.0008 \ 
-.001 1 

-.002 
-4.55 

-4.66 / 



We have summed over all n such that kp.z + nQ lies within 
the first Brillouin zone of the superconductor (|n| < 8). Ma- 
trix elements evaluated at the Fermi momenta of other filled 
bands will be similar; there will be small differences when 
kp.j is such that other n are allowed, but these will ultimately 
result in a very similar induced gap. We do not need to worry 
about k x ^ k' x because if k' x = k x + Q, then only one of 
them can be near the Fermi surface of the nanowire. Since 
the matrices are effectively diagonal, we can calculate the in- 
duced gap from Eq. ( |25] >, shown in Fig [7] as a function of the 
tunneling strength t, along with the induced gap for a wider 
wire with w = lOOnm. 
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